THE CONE PERCOLATION ON 
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Abstract. We study a rumour model from a percolation theory 
and branching process point of view. The existence of a giant 
component is related to the event where the rumour spreads out 
trough an infinite number of individuals. We present sharp lower 
and upper bounds for the probability of that event, according to 
the distribution of the random variables that defines the radius of 
influence of each individual. 



1. Introduction 

We study long range dependent oriented percolation processes on a 
tree through its most basic propriety: the existence of a giant connected 
component. The starting point for approaches to rigorous percolation 
theory beyond the nearest neighbor independent setup on is due 
to several authors around the nineties. Grimmett and Newman [5] in 
1990 study percolation on x Z, Burton and Meester [3j in 1993 study 
phase transition for a long range independent percolation model on a 
stationary point process in M'^, Lyons [S] put out the first version of his 
book Probability on Trees in 1994 while Benjamini and Schram [1] in 
1996 have they Percolation beyond Z'^ published, just to name a few. 
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Lebensztayn and Rodriguez [7] in 2008, propose a model on general 
graphs named disk percolation where a reaction chain starting from 
the origin of the graph, based on independent copies of a geometric 
random variables with parameter q G [0,1], defines the existence or 
not of a giant component. They obtain a sufficient condition for the 
existence of phase transition based on q, which means the existence of a 
non-empty subcritical (no giant components) and supercritical (giant 
components with positive probability) phases. They associate their 
model to a rumour or an epidemic process. In this paper, instead of 
working in a general family of graphs we focus on homogeneous trees 
and instead of fixing the random variable which defines the radius of 
infection or the radius of infiuence of each vertex to be geometric, we 
consider general random variables. So, as a result, instead of having a 
phase transition phenomena dependending on a point in a parametric 
space, we have that phenomena depending on the family of general 
positive random variables. 

We consider a process which allows us to associate the activation 
dynamic on the set of vertices to a discrete rumour process. Indi- 
viduals become spreaders as soon as they heard about the rumour. 
Next time, they propagate the rumour within their radius of influence 
and immediately become stifiers. Our main interest is to establish 
whether the process has positive probability of involving an infinite 
set of individual. Besides, we present sharp lower and upper bounds 
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for the probability of that event, according to the general distribution 
of the random variables that defines the radius of influence of each 
individual. We say that the process survives if the amount of vertices 
involved is infinite. Otherwise we say the process dies out. 

Consider the homogeneous tree such that each vertex has d-\-l 
neighbours, d > 2. Let V{Tfi) the set of vertices of T^. We single out 
one vertex from V(Td) and call this O, the origin. For each two vertices 
u,v E V(Trf), we say that u < v if u belongs to the path connecting O 
to V. Besides, for two vertices u,v such that u < v let d{u,v), be the 
distance between u and v, as the number of edges the path from u to 
V has. Now, let us define 

T+(u) ^{ve V{Td) -.uKv}. 

Pick a ^; e V{Td) such that d{0,v) = 1 and consider T+ = Td\T+(w). 
For Tj" we define 

dT+{u,n) = G T+ : d{u,v) = n}. 

Now we define the Cone Percolation Process in T^. Let {Rv}{vev{Ta}} 
and i? be a set of independent and identically distributed random vari- 
ables. We define Pk = P[R = k] for k = 0,1, . . . To avoid trivialities we 
assume po e (0, 1). Besides, for each u e V{Td), we define the random 
sets 

B^^{ve V{Td) : d{u, v) < R^} (LI) 
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and consider the non-decreasing sequence of random sets /q C /i C • • • 
defined as Iq = {O} and inductively J„+i = IJug/„ for all n > 0. 
Let / = IJn>o be the connected component of the origin. Under the 
rumour process interpretation, / is the set of vertices which heard about 
the rumour. We say that the process survives if |/| = oo, referring to 
the surviving event as V. 

Consider P+ and P the probability measures associated to the pro- 
cesses on and (we do not mention the random variable R unless 
absolutely necessary). By a coupling argument one can see that for a 
fixed R 

p+[i/] < p[y] (1.2) 

By the other side, by the definition of and its relation with 
we have that for a fixed R 

P+[y] = if and only if F[V] = 0. (1.3) 

The paper is organized as follows. Section [2] presents the main re- 
sults. Section |3] brings the proofs for the main results together with 
auxiliary lemmas and handy inequalities. Section H] presents results for 
the heterogeneous setup of the Cone Percolation Process. Finally, in 
Section [5] we present examples where some conditions can be verified. 
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2. Main Results 

Theorem 2.1. Consider the Cone Percolation Process on with 
radius of influence R 

(I) IfE{d^) >l+po then, P+[1/] > 0, 
(II) IfE{d^) < 2 - i then, = 0. 

Let p and ip be, respectively, the smallest non-negative root of the 
equations 

E(p<^") + (1 - p)po = p, (2.1) 
E(V'd^('^''-i)) = i;. (2.2) 

Theorem 2.2. Consider the Cone Percolation Process on Tj. Then, 

1- p< P+(y) < 1 - -0- 

Theorem 2.3. For the Cone Percolation Process on Td with radius of 
influence R, it holds that 

1- (l - p'-^) Po-E (p^'^") < F[V] < 1-E (V;W(^"-^)) . (2.3) 

3. Proofs 

3.1. Auxiliary Processes. Let us define two auxiliary branching pro- 
cess, being the first one {Xn}n€f^- For this process, the associated ran- 
dom variable is X, assuming values in {0,d,d^, . . .} such that 

P[X = 0]=p„, 

PlX^d''] for A; = 1,2,... 
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whose expectation is 

E[X] = E[rf«]-po (3.1) 

and whose generating function is 

Ms) = E[s^] = E[s'''] + (1 - s)po. (3.2) 

The second auxihary process is {yn}nm- For this process, the associ- 
ated random variable is Y, assuming values in {0, d, d+d'^, . . . , Yl'i=i 
such that 

p\y = ^^^^] =PkioTk = o,i,2,... 

I ct — 1 J 

whose expectation is 

E[Y] = ^(E[rf^] - 1) (3.3) 
and whose generating function is 

Ms) = E[s^] = E[s^(''''-^)]. (3.4) 

3.2. Proofs. 

Proof of Theorem \2.1\ 

By a coupling argument one can see that our process dominates {Xn}nm- 
This process survives as long as E[X] > 1 therefore from (13.1 p our pro- 
cess survives if E[c?-'^] > 1 + po; proving (I). 

By the other side, also by a coupling argument, our process is domi- 
nated by {3^„}„gN. That process dies out provided E[y] < 1 therefore 
from f l3.3p our process dies out if E[(i^] < 2 — i, proving (II). □ 
Proof of Theorem \2.^ 
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In order to find the extinction probability of {Xn}nen (Grimmett and 
Stirzaker( [HI p. 173]), let us consider the smallest non- negative root of 
the equation p = Lpx{p)- Therefore from (I3.2p 

E[/"] + (l-p)po = P 
and by construction of the processes, as P+[y^] < p, we have that 

l-p<P+m. 

In order to find the extinction probability of {3^n}nGN (Grimmett and 
Stirzaker( [6l p. 173]), let us consider the smallest non-negative root of 
the equation -0 = lpy{'<P)- Therefore from (13 ■4p 

E[^^('^"-i)] = ^ 
and by the construction of the processes, as P+[y^] > -0, we have that 

p+[y] < 1 - ^/^. 

□ 

Proof of Theorem \2.3\ 

Observe that except for the root, all vertices see towards infinity a 
tree like T^. So, assuming Rq = k the probability for the process 
to survive is larger or equal than the probability of the process to 
survive from at least one of the d^~^{d+ 1) trees that have as root the 
furthest infected vertices. By the other side, still assuming Rq = k, 
the probability for the process to survive in T^^ is smaller or equal 
than the probability for the process to survive from at least one of the 
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{d + l){d^ — l){d — 1)^^ trees like that are seen from each active 
vertices by its own, independently. So, ior k = 1,2, . . . 

1 - (1 - P+[y])('^+i)'^''' < ¥[V\Ro = k]<l-{l- ¥+[V])^^'^'-'\ 

From this and from Theorem 12.21 follows (12. 3p . □ 

4. Heterogeneous Cone Percolation on T| 

Suppose we have two sets of independent random variables, {Rz}{zef>i} 
and {Rv}{y^v(j+)}, such that for all 2; G N and all m G V such that 
d{0, u) = z, Ru and Rz are equally distributed. Besides assume 
F[Rz = 0] < 1 for all zeN. 

We define the Heterogeneous Cone Percolation Process from the set 
of Bu defined in (11. ip . For n G N fixed and u < v G V(T^), consider 
the event 

V^^ : Process starting from u reaches v in at most n steps. 

For a fixed integer n, let = {O}. Besides, for j = 1,2, .. . we 
define 

X;= U {vedT^-.V:^, occurs}. 
Again, for all j = 1, 2, . . . consider 

z; = |x;|. 

So, for all fixed positive integer n, {^"}j>o is a branching process 
dominated by the number of vertices v G dTj^ which are activated. 
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Lemma 4.1. Consider n fixed. For nj := Ei[Zj-], the mean number of 
offspring on generation j for the process {Zj'}j>o, it holds that 

where = P[T4"^], for any fixed pair u < v such that d{0,u) = jn 
and d{0, v) = {j + l)n. 



Proof of Lemma 4 ■ 1 



For fixed j and n, consider dT^ = {ui,U2, ■ ■ ■ ,Udr^}. So we can write 
Z" as Yli=i -^{v;"„,}- Taking expectation in both sides finishes the proof. 

Lemma 4.2. Consider n fixed and = P[V^^], for any fixed pair 
u < V such that d{0, u) = jn and d{0, v) = {j + l)n. 

n— 1 k 

p] > -l[m,k+^ < A; + 1 - i]]. 

k=0 i=0 



Proof of Lemma \4-S\ 

For any fixed pair u < v such that d{0, u) = jn and d{0, v) = [j + l)n 
we have that 



and so 



n-l 

n 

k=0 



[j{Rjn+i >k+l-l} 



i=0 



p n 



\}{R3n^^ > A; + 1 - i} 



n-l 



.i=0 
k 



> l[^[[j{RJn+^>k+l-^}]. 



k=0 \i=Q 
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The inequality is a consequence of the FKG inequahty (Grimmett [U 

p.34]). □ 

Theorem 4.3. The Heterogeneous Cone Percolation Process in has 
a giant component with positive probability if for some fixed n, 

n— 1 k 

hminf d'^TT[l - TTP[i?,fc+i <k+l-i]]>l. (4.1) 



1— >oo 

fc=0 i=0 



Proof of Theorem \4-^ 

From Theorem 1 in Souza & Biggins ([21 p.39]), a branching process in 
varying environments is uniformly supercritical if 

lim inf /ij > 1 . 

From Lemma 14.11 and Lemma 14.21 that is what happens if (14.11) holds. 
From the fact that 

1 

^ < — 



one can see that the Heterogeneous Cone Percolation Process has a 
giant component with positive probability if 

n— 1 k 

lim inf TT[1 -W¥[Rjk+i <k + l-i]]>l. 

fc=0 j=0 

□ 

5. Examples 

Example 5.1. Consider a Cone Percolation Process in T^^, assuming 
V[R = l]=p= 1-P[i? = 0]. 

In words R ~ B{p). 
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• If p > d-^ then, F[V] > 0. 

• If p < then, F[V] = 0. 

By the definition, one can see that 

P[F^] = (i_p) + (p^[i/-])'^+ip. 

Observing that the upper and lower process presented by {Xn}nm 
and {yn}nm presented in session [3TT] are the same, we see that 

being ip the solution of 

pijj'^ - ■?/' + l- p = 0. 

Example 5.2. Consider a Cone Percolation Process in Td, assuming 
P(i? = ^) = (1 _ p^pk^ A: = 0, 1, 2, . . . 

In other words R ^ Q{1 — p). From Theorem 12.31 

• If dp"^ -2dp+l <0 then, F[V] > 0. 

• If 2pd < 1 then, F[V] = 0. 

As a consequence of this and f ll.Sp . for fixed 



^<inf{p:P[V^]>0}<l-Y/l-^. 



Example 5.3. Consider a Cone Percolation Process in T^, assuming 



p(/2 = fc) = ^ = o,i,...,n 

If + 1 - p)" - (1 - p)" > 1 then, F[V] > 0. 
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• li2d- d{pd + 1 - p)" > 1 then, ¥[V] = 0. 
Let = 2 and R ~ i3(4,i). 
Therefore p and i/j are, respectively, solutions of 

+ 4a;8 + Qx^ + 42;2 _ 16x + 1 = 0, 

+ 4x^^ + + 4a;2 - IGx + 1 = 0. 

So p = 0.0635146 and ip = 0.06350850, which implies that 

0.937435919 < F[V] < 0.937435962. 

Let d = 4andi?~S(4, 1). 

Therefore p and i/j are, respectively, solutions of 

x^^^ + 12x^^ + 54x^^ + lOSx^ - 256x + 81 = 0, 

+ 12x^^ + 54x2° + lOSx^ - 256a: + 81 = 0. 

So p = 0.3208787235 and ^ = 0.3208787200, which implies that 

0.682158629 < F[V] < 0,682158630. 

Example 5.4. Consider a Heterogeneous Cone Percolation Process on 
T^, assuming that Rj are Bernoullis, that is, 

F[Rj = 1] = 1 - F[Rj = 0] for j = 0, 1, 2 . . . 

By applying Theorem 14.31 with n = 1 one can see that the Heteroge- 
neous Cone Percolation Process on survives with positive probabil- 
ity if 

liminf (iP[/?j = 1] > 1. 
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